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The scaling behavior of cyclical growth {e.g. cycles of alternating deposition and desorption 
primary processes) is investigated theoretically and probed experimentally. The scaling approach 
to kinetic roughening is generalized to cyclical processes by substituting the time by the number of 
cycles n. The roughness is predicted to grow as where jS is the cyclical growth exponent. The 
roughness saturates to a value which scales with the system size L as L°' , where a is the cyclical 
roughness exponent. The relations between the cyclical exponents and the corresponding exponents 
of the primary processes are studied. Exact relations are found for cycles composed of primary linear 
processes. An approximate renormalization group approach is introduced to analyze non-linear 
effects in the primary processes. The analytical results are backed by extensive numerical simulations 
of different pairs of primary processes, both linear and non-linear. Experimentally, silver surfaces are 
grown by a cyclical process composed of electrodeposition followed by 50% electrodissolution. The 
roughness is found to increase as a power-law of n, consistent with the scaling behavior anticipated 
theoretically. Potential applications of cyclical scaling include accelerated testing of rechargeable 
batteries, and improved chemotherapeutic treatment of cancerous tumors. 
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I. INTRODUCTION 



From DLA to MBE, kinetic models of growth and aggregation have attracted much attention Q| [D @ (|] in 
the last two decades due to diverse interests in physics, biology, chemistry and engineering. Kinetic roughening of 
non-equilibrium surface growth were of particular interest. Crystal growth, the growth of bacterial colonies, and 
the formation of clouds in the upper atmosphere |^ are all examples of non-equilibrium phenomena which grow 
self-affine rough surfaces. On a fundamental level, surface growth problem is a paradigm for a class of problems 
in non-equilibrium statistical mechanics. One crucial aspect of this class is the signature of scale invariance and 
universality, very similar to those observed in the equilibrium critical phenomena, or in non-linear dynamical systems 
. Investigations focused primarily on the scaling behavior of the surface roughness , , |^ (|] • More recently 
they have touched upon other aspects such as the distribution of the surface width , distributions of the height 

and the average height velocity |Q| , density of extrema |l2) , persistence and the maximal height |0 . 

Processes which generates rough surfaces can be divided into two classes: (i) growth processes {e.g. by deposition, 
absorption) where the surface grows by adding material; and (ii) recession processes {e.g. by erosion, dissolution, or 
desorption) where material is being taken off such that a rough surface is generated. Examples of class (i) include 
crystal growth Q ||^, electro-platin g [|l5[ , and biological growth j|] |^ which have been studied widely. 

Class (ii) includes chemical dissolution |17||, and received much less attention than the first class. 

Self-affine surfaces generated by growth (or recession) can be described using scaling analysis of the surface rough- 
ness. In general, the surface is characterized by the height h{f, t) appropriate to the {d = d—1) dimensional substrate 
of linear size L. The width of the surface W{L, t) at a time t is given by 

W{L,t) = ([h{r,t)-{h{r,t))yy, (1) 

where {h{r,t)) is the average height which is linear in time {h{r,t)) = vt, v being the average growth (or recession) 
velocity. The angular bra-ket () denotes average over both lateral sites and the ensemble of surface configurations. 
W{L,t) scales as 

W{L, t) ^ f{Lim) - f{L/t'/n, (2) 

where ^(t) ~ t'^^'^ is the lateral correlation length and / is the scaling function. 
For large time t:^ L"": W^L", 
while for t<i^: W ^ tl^ , 

where f3 — a/ z is the growth exponent. The height difference correlation function 

A(r,t) = (^(^h{f^ + f^t)-h{fo,t)y), (3) 
obeys scaling similar to that of the roughness 

A(r, t) ^ r'V{r/m) ^ r^V{r/t'^n- (4) 

For r<ii/^: A ~ r^", 

while for r>ii/^: A - t^/S^ 

The relation between the height-difference correlation function A(r, <) defined above and the equal-time height-height 
correlation function C(r, t) — (h{fo + r, t) h{ro, t)) is simple 

A(r,t) = 2[C(0,t)-C7(r,t)]. (5) 

Dynamic Scaling in Fourier space is easily studied in terms of the structure factor S{q,t) — (^h{q,t) h{~q,t)'j, where 
h{q,t) is the Fourier transform of the height h{f,t). The scaling hypothesis Eq. (2) can be translated to the Fourier 
space such that 

S{q,t)=q-''-''^g{q/t-^^% (6) 

Surface width W can be readily calculated from S{q,t) using the relation W'^{L,t) = {^/L'^)J2q^{'lT^)- 
theoretical analysis of surface growth it is convenient to work in Fourier space and compute first the structure factor 
rather than the roughness W itself. Also, in some experiments the surface is probed by scattering processes which 
provide its structure factor. 
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All processes within the same universality class share the same critical exponents. Their continuum growth equations 
differ at most by terms which are rendered irrelevant by the renormalization group flow. The asymptotic continuum 
stochastic equations corresponding to different universality classes (indexed by i = 1,2...) are of the form of a Langevin 
equation 

^^=A.m + ,y,(r,t)+x.„ (7) 

where Ai{h} is a local functional depending on the spatial derivatives of h{r,t) and the noise r]i{r,t) reflects the 
the random fluctuations in the deposition process and satisfies {r]i{'r,t)) = and ti) 77^(7^2, ^2)) = 2DiS'^{fi — 

^2) S{ti — ^2)- Specific generic processes and their universality classes are reviewed later (Sec. III). 

In the present paper we focus on cyclical processes in which deposition and desorption are occurring alternatingly. 
Examples of cyclical processes are abundant in nature and technological applications. A technologically important 
example are rechargeable batteries, where metal is electrodeposited on an electrode during the discharge, followed 
by partial electrochemical dissolution of this metal during recharging. In chemotherapeutic treatment of cancer the 
malignant cells are subjected to a recessive cyclical process. 

The basic premise of our scaling approach to cyclical processes, is that the number of cycles n should substitute for 
the time variable. So we propose for cyclical processes the scaling relation of Eq. (1) be replaced with 

WciL, n) ^ UL/Ut)) - UL/n^'^l (8) 

where ^ "n^^^ is the correlation length and /c is the cyclical scaling function. 
For large n 71 » L^: W ~ L", 

while for n <^ L^: W ^ , 

where f3 = a/z is the growth exponent. 

The rest of the paper is organized as follows: The next section II is devoted to our analytical analysis. We will 
prove that the equation (8) for cyclical scaling holds asymptotically for linear processes. In addition we will show 
how to obtain the scaling exponents of the cyclic process, given those of the two primary processes. An approximate 
RG approach is introduced to study non-linear effects. Sec. Ill contains a brief review of the generic universality 
classes and their lattice algorithms. Then we introduce new algorithms to implement numerically different recession 
processes. Actual results of our simulation of cyclical processes are described in Sec. IV. In Sec. V the experimental 
findings from cyclical electrodeposition/dissolution of silver are discussed. Sec. VI contains the conclusions with a 
discussion of potential practical applications of the cyclical scaling approach. A short summary of some of the results 
was published in Ref. 



II. ANALYTICAL RESULTS 



The analytical approaches are all based on the stochastic equation, like Eq. (7), which describes the growth process. 
Thus, we begin by obtaining the stochastic equation of the cyclical process. The first primary process is denoted by 
i — 1 and the second hy i — 2. The durations of the first and the second processes are Ti — pT and T2 = (1 — p)T, 
respectively (p and 1-p are the fractional parts of both processes). The total time period for one cycle is T = Ti +T2. 
The cyclic growth equation in terms of the basic two processes can be expressed as 

dh 

= [a^h + 771 + v^Mp - fit)) + [a^h + 7?2 + V2Mf{t) - p), (9) 

ot 

where f{t) is defined as the fractional part of t/T and Q(x) is the unit step function. 



A. Linear Primary Processes - Exact Calculation of Sc{q,t) 

First we will consider cyclical processes for which both the primary processes are linear. For linear processes 
Ai{h} — ai{'S7)h(r,t), where ai(V) is a linear differential operator. The time reversal symmetry is obeyed in this case 
if the height is measured relative to the average height. The average height indeed depends on the growth {vi > 0) 
and the recession (vi < 0) nature of the primary processes. In terms of the average velocity Vc = pvi + (1 — p)v2, it is 
given by 

{h{f, t))/T = nv, + f{t)Q{p - fit)) + [{v, - V2)p + 1^2 (/(<))] e(/(i) - p). (10) 
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The roughness is insensitive to the sign of Vi and hence will not distinguish between a growth/growth and a 
growth/recession cyclical processes as long as the basic processes remain linear. 

For linear processes Langevin equations of the form of Eq. (7) are easily solved in Fourier space to yield (assuming 
spatial isotropy in the basal plane) the structure factor 

S{q,t) = exp{-2a(9)t}5(9,0) + ^[1 - exp(-2a(g)0], (11) 

where S{q, 0) is the structure factor at f = which contains the information of the initial roughness. During the 

nth cycle of the cyclical growth, the structure factor Sc{q, {{n — 1) + p)T^ generated by the first primary process (of 
durationTi = pT) is assigned as the initial condition for the second primary process. The second process lasts for 
T2 = (1 — p)T to yield the structure factor Sc{q, nT) of the cyclical process after n cycles. This is again used as the 
initial structure factor for the first process in the (n + \)th cycle. 

During the first cycle, the structure factor after the completion of the first primary process becomes 

5(g,ri) =exp{-2ai}5(g,0) + -^[l-exp(-2ai)], (12) 

where di = ai{q)Ti was defined. This is the initial structure factor for the second primary process and after the first 
complete cycle we obtain 



S{q, T2) = exp{-2(ai + 02} S{q, 0) 

+ exp{-2a2}| — [1 -exp(-2ai)]| + — [1 - exp(-2a2)], (13) 

L fll J 012 



where Sc{q, T) = S{q, T2). Proceeding in this manner, we finally arrive at the structure factor Sc{q, n) = S{q, nT) 
of the cyclic growth after n cycles as a geometric series which can be readily summed to yield 



Sc{q,n) = exp{-2acn}S{q,0) 

1 — exp(— 2acn) 



+ 

where = UrT with 



— - exp(— 2a2){l — exp(— 2ai)} H ^{1 — exp(— 202)} 

fll 02 



1 — exp(— 2ac) 



(14) 



ac = [aip + a2{l - p)]. (15) 
In the scaling limit of small q, Eq. 14 for Sc{q,n) reduces to 

Sc{q, n) ~ [1 - exp ( - 2ae(g) Tn)], (16) 

ctcyq) 

with the effective noise strength for the cyclic process is defined as 

D,=pD, + {l-p)D2. (17) 

In terms of the effective parameters ac and Dc, the above structure factor for the cyclical process resembles that 
of a generic linear growth process (see Eq. (11) with the number of cycles n as the new time variable. Effectively the 
time is being coarse-grained over a period of T by eliminating the high frequency (> 2tt/T) modes. Hence we can 
use the standard scaling analysis (Eq. (6)) of the structure factor to determine the scaling exponents in the case of a 
cyclical growth. 

After a large number of cycles n the interface width becomes saturated. In that limit, Eq. (15) becomes Sc{q, n) ~ 
■ The roughness exponent of the cyclic process is determined by the q Q divergence of . Since (g) ~ q^^ , 
it is the process with smaller Zi which dominates the asymptotic cyclical roughness, and the roughness exponent is 
given by ac = min{a\,a2) = \{min{z\, Z2) — {d— 1)}. The primary process with the smaller roughness imposes its 
roughness exponent on the combined cyclical process. The larger ai appears as a correction to the scaling exponent. 
Whether or not it is the leading one depends on how its contribution compares with that of the subleading term 
in ai{q) of the dominating primary process. Note that a subleading term in ac{q) might affect the behavior on a 
smaller scale, if its amplitude is large. In that case, the leading behavior takes over only beyond a crossover length 
(at which both contributions are comparable). Since the amplitudes of ai and a2 in Sc{q,n) are proportional to p 
and (1 — p), respectively, the longer the non-dominant process lasts, the larger is the crossover scale, as could be 
expected. However, although in ac{q) only the dominant ai{q) is important beyond this crossover scale, this is not 
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the case for the effective noise-correlator Dc, which is a scale independent constant in Eq. (16). Thus, the amplitude 
of the leading power-law roughness is determined by both the primary processes. 

In the growing phase of the interface roughness, the dynamic exponent z will dictate the cyclical power law behavior. 
Since n is multiplied by ac = [aip + 02(1 — p)\T in Eq. 16 and ai{q) ~ the process with smaller Zi will again 
dominate the cyclical dynamics in the q ^ Q asymptotic limit to yield Zc = min{z\,Z2)- Then for the initial cycles 
(nT <C L^), the surface width will scale as , with fic = ctcl Zc- In essence, the scaling exponents of the cyclic growth 
will be identical to those of the primary process with the smaller Zi. 



B. Coarse-Graining Approach: The CycHc Propagator Gc{r,t) 

In the previous section the structure factor Sc{q,t) (and hence its Fourier Transform (FT) Sc{r,t)) was obtained 
exactly. The same method of successive integration of the cyclical equation of motion may be applied to derive the 
cyclic propagator Gc{r,t) (or its FT Gc{q,t)). However, the expression is cumbersome and not very useful. Since we 
only need the long-time behavior, we will take a different route of coarse-graining the equation of motion such that 
the equation for time scales larger than T will be derived and Gc{r, t) could be read from it. This will also provide a 
direct connection to the RG approach introduced in the next section to study the behavior of non- linear systems. 

Our starting point is Eq. (9), which we choose to integrate over one cycle, such that the remaining equation will 
be a difference equation between the average heights of the consecutive cycles. So assume we average Eq. (9) on time 
t € [nT, (n -I- 1)T]. The integral will be divided into integration over the two intervals: interval (1) [nT, nT + Ti] and 
interval (2) [nT + Ti,{n + 1)T]. We define the average height in the n*'' cycle as 



1 



n+1 



H{r,n) = - dth{r,t). (18) 

^ J n 

Our goal is to obtain the equation of motion for i/(r, n) on time scales t > T. The FT of Eq. (9) is 

= [ai{q) h{q, t)+m+ v^]e{p - /{t)) + [a2{q) h{q, t)+V2 + V2Mf{t) - p) 
= t) + [Aai(9) e(p - f{t)) + Aa2{q) e{f{t) - p)] h{q, t) 

+{Vi + vi) e{p - fit)) + (r?2 + ^^2) e(/(i - p)), (19) 

where Aoi(g) = ^ai(g) — ac{q)^ and Aa2(g) = ^02(9) — adqij- Let us now integrate over one cycle each of the 
terms, beginning with the l.h.s. 



-T r 



dt = ^ [h{q, (n + 1)T) - h{q, nT)] 

= ^[H{q,n+l)-H{q,n)] 



+ ^ { (h{q, (n + l)r) - H{q, n+1)) - (h{q, nT) - H{q, n)) }. (20) 

The first term can be expressed as ^j/^"'* ■ The second term contains the differences between the height at the 
beginning of the cycle and its average over a cycle. This term thus reflects the behavior within each cycle and is thus 
irrelevant to the behavior on the coarse grained scale and will be dropped. (Note that it will vanish asymptotically 
and it is irrelevant since in the continuum limit it features a second derivative with respect to t). 

The integration of the first term on the right-hand side yields 



n+1 

aM h{q, t) dt = ac{q)H{q, n). (21) 



Integrating the second term 



1 

T 



n+1 

[Aai(g-) e{p - fit)) + Aa2{q) e(/(t) - p)] h{q, t) dt 



h{q,t)dt+^^ h{q,t)dt 
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^Aai{q)p[—J h{q,t)dtj+Aa2{q)il-p)[—J hiq,t)dtj 

= Aai {q) p Hi {q, n) + Aa2 (q) {1 - p) H2 {q, n) , 
= [Aai(g> + Aa2(g)(l - p)]H{q, n) 

+{AaMp[Hi{q,n) - H{q,n)] + Aa2{m - P)[H2{q,n) - H{q,n)]}, (22) 

where Hi{q,n) and H2{q,n) are the average heights during the two primary processes respectively (on the n*'' 
cycle). Upon coarse-graining 



n+l 

[Aaiiq) e(p - fit)) + Aa2{q) Q{f{t) - p)] h{q, t) dt 
= [Aai(q> + Aa2(g)(l-p)]i?(<f,n) (23) 



This term contains information on the structure factor within the cycle and must vanish in the coarse-grained 
equation, since ac was chosen (Eq. (15)) such that 

(ai-ac)j3+(a2-ae)(l-p)=0. (24) 

The velocity term yields the average velocity 



dt{vi Q{p - fit)) + V2 e(/(i) - p)) = V,. (25) 



The noise 77's are random variables and we define its coarse-grained term by 



1 



n+p rn-\-l 



ryc(n) = -=[/ dtrn{t)+ dtTj2{t)]. (26) 



n+p 



It obeys {ridn)) = and {ric{n) rjc{m)) = 2DcSn,m, and Dc satisfies Eq. (17). 
Collecting all the terms, the coarse-grained difference equation for H{q, n) is 

= aMm, n) + rjM ") + Vc (27) 

1 An 

For 1, the difference equation is equivalent to the corresponding differential equation ~ 
The cyclical propagator may be read from this linear equation 

Gc{q, n — m) = exp{— ac((?)(n — m)}Q{n — m). (28) 

This propagator obeys the usual relation to the structure factor if n is treated as a continuous variable replacing 
the time t 

SMn)^^[l-\GMn)?]. (29) 



C. Non-Linear Primary Processes 

Stochastic non- linear equations of that type can be analyzed using the perturbative dynamic RG approach |^ [ p3| . 
We develop an approximate RG technique for cyclical processes with one or both the primary processes being non- 
linear. The first step is to set aside (for the initial RG iterations) all the non-linearities from the primary processes 
and take only the linear Langevin equations. It was shown above how to solve for the structure factor of such cyclical 
processes and get an effective S'c('Z, n) (see Eq. (16)) of a non-cyclic process. Similarly we have shown in the previous 
subsection how to obtain the propagator Gc{q,n) of the same effective linear, non-cyclical process. The second step 
is to take the effective "free" propagator Gc{q, n) and then add back all the non-linear terms of the primary processes 
as perturbations. The bare couplings of the non-linear terms are multiplied by p and (1 — p) to take into account the 
relative durations of the two primary processes. The third and the final step is to study the RG flow of the parameters 
and determine the fixed points of the transformations and hence the scaling exponents. This dynamic RG procedure 
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is approximate in the sense that the initial flow of the couplings is shifted, but as long as the starting point is not 
close to a separatrix in the parameter space {i.e. a border line between basins of attraction of two different fixed 
points) this will not alter the ultimate fixed-point of the RG flow. 

The cyclical process might have only one non-linear primary process with the non- linear perturbation being relevant 
under RG transformations. Then the cyclical scaling exponents are given by the exponents of the non-linear primary 
process. Otherwise, the non-linear term turns out to be irrelevant with respect to the dominant linear term (present in 
the coarse grained cyclic free propagator) of the other (linear) primary process, which force its own scaling exponents 
for the cyclic growth. 

Both the primary processes might contain non-linear terms which are treated as perturbations to the free propagator 
as described above. In most cases, one of the non-linear term will make the other one irrelevant under the RG flow 
and the primary process containing the relevant non-linearity will carry through its scaling exponents for the cyclical 
process. If both the non-linear terms are relevant, then the possibility of new cyclical growth exponents (different than 
those of both primary processes) can not be ruled out. A prime candidate for such a behavior will be a cyclical process 
in d = 2 + 1 in which one process has a relevant non-linearity while the second process has a relevant anisotropy in 
the basal plane 

III. GROWTH MODELS AND UNIVERSALITY CLASSES 

In this section We present a brief review of the basic models of kinetic growth. The discrete lattice models for 
simulations corresponding to each universality class are presented. We were able to generalize most of the lattice 
models for desorption processes (reverse of growth) to be used in absorption/desorption cycles. These new desorption 
algorithms are described in detail. 

1) Random Deposition (RD) - This is simplest of all possible growth processes. From a randomly chosen site 
from the surface, a particle falls vertically until it reaches the top of the column under it, whereupon it is deposited. 
In this case 

Ard = 0. 

Scaling exponents are: /? — 0.5 in all dimensions but a is not defined for this model because the interface is never 
saturated. 

RD Algorithm - A column i is chosen (d=l-|-l) randomly and its height h(i,t) is increased one. For the reverse 
process we just decrease the height h(i,t) of an arbitrarily chosen site i by one. 

2) Edwards — Wilkinson (EW) Universality - In this model we have surface relaxation (in addition to the 
random deposition) which is introduced by the term 

Aew = i^V^/i. 

We can exactly calculate the scaling exponents: a = /3 = and z = 2. In EW growth a randomly 

deposited particle can diffuse along the surface up to a a finite distance and sticks to a local height minimum. Due to 
this relaxation, surface becomes smooth compared to the random growth model and finally the interface roughness is 
saturated because of the correlations among the neighboring heights, to a value ~ L", where L is the lateral size. 

Family Model - This model was introduced by Family |2^ to simulate EW growth. A particle is dropped on a 
randomly chosen column i (in d = 1 + I) and sticks to the top of the column i, i-l-1 or i-1, depending on which of the 
three columns has the smallest height. To simulate a desorption process (with EW exponents), the height h{i,t) of 
the site i is compared to the heights h{i— l,t) and h{i + l,t) of the neighboring sites. Then we simply take a particle 
off from the column i, i-|-l or i-1, depending on which of the three columns has the largest height. If case of a tie 
involving the site i we take out a particle from that site, otherwise the tie is broken randomly with equal probability. 
The whole process can be thought as desorption with surface relaxation. The relaxation length is restricted to the 
nearest neighbors because the scaling exponents are independent of the relaxation length. In Fig. 1 the roughness 
W is plotted against time on a log scale for increasing system sizes. The slope of the straight line fitted to the early 
time roughness yields the growth exponent /? = 0.23 ± 0.02. The roughness exponent a = 0.48 ± 0.03 is extracted 
from the dependence of the saturation roughness on the linear size of the system (Inset of Fig. 1). These values 
of the scaling exponents are in agreement with the corresponding EW values. Indeed, the surface roughness for the 
desorption process shows (Fig. 1) very similar behavior to that obtained for the growth process So the discrete 
rules described above can be treated as a valid algorithm for EW desorption. 

3) Kardar — Parisi — Zhang (KPZ) Universality - This describes growth in a direction locally normal to the 
interface |2^. Its leading effect is to add a non- linear term to the EW surface relaxation term 
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Scaling exponents for KPZ equatfon are exactly known in d = 1 + 1: ct = 1/2, (3 = 1/3, and z = 3/2. In d = 2 + 1 
approximate (from numerical simulations) values for the exponents are: a ~ 0.39, f3 ~ 0.24 and z ~ 1.61. KPZ can 
be simulated using different atomistic growth algorithms of which two we will describe because we generalized those 
to the case of desorption and will be used in our simulations of cyclical growth. 

Ballistic Deposition (BD) Model - A particle is released from a randomly chosen position above the surface, located 
at a distance larger than the maximum height of the interface. The particle follows a straight vertical path until it 
reaches the surface, whereupon it sticks ||2^. If h(i,t) is the height of the column i (chosen randomly in d = 1 = 1) 
at time t then the BD growth rule is: h(i,t+l)=max[h(i-l,t), h(i,t)+l, h(i+l,t)]. For the reverse process the algorithm 
will be changed to h(i,t+l)=min[h(i-l,t), h(i,t)-l, h(i+l,t)]. Although physically unrealistic (contrary to its growth 
counterpart) this desorption model is formally the "anti-BD" process. 

Restricted-Solid-On-Solid (RSOS) Model - This algorithm (also known as KK model) introduced by Kim and 
Kosterlitz gives KPZ exponents and known to yield reliable results even for small system sizes. The growth 
rule is to randomly select a site on a cubic (d-1) dimensional lattice and to permit growth by letting the height of 
the interface hi — > /i^ + 1 provided the restricted sohd-on-solid restriction on neighboring heights \Ah\ = 0,1,...., N 
{N > 1) is obeyed at all stages. In a similar way we can simulate an erosion process where we decrease the height 
{hi ^ hi — 1) of the site i provided the RSOS condition |A/i| = 0, 1, N is satisfied. In Fig. 2 we show the results 
of simulations using this desorption rule. The values of the scaling exponents {(3 = 0.33 ± 0.02 and a = 0.52 ± 0.02) 
obtained are consistent with the corresponding KPZ values. 

The RSOS deposition or desorption rule described above is defined starting from a flat interface at t = 0. In a 
cyclical surface growth model with a primary process not obeying the RSOS condition can destroy the height difference 
restriction of the RSOS model. We can still use a growth rule similar to the RSOS model described above, which 
seems to be behaving like a KPZ growth (as far as the scaling exponents are concerned). In this extended model 
of RSOS, we choose a site i randomly and add a particle on it only if the height of that site is less than or equal 
to the heights of the neighboring sites (which would give us the normal RSOS with \Ah\ = 0, 1 starting from a flat 
interface). For the reverse process a particle is taken off from a site only if the height of that site is greater than or 
equal to the heights of the neighboring sites. 

4) Mullins — Herring (MH) Universality - In conserved growth situations where " surface diffusion" is dynami- 
cally significant in the absence of any EW relaxation process, the growth process may belong to the MH universality 
class (also known as Das Sarma-Tamborenea (DT) or Wolf- Villain (WV) class) |^^, j29j. The linear surface 
diffusion equation for MH has a term 

Amh = -KV^h. 

The critical exponents for the MH growth universality are exactly known theoretically: a = (5 — d)/2, (3 = {5 — d)/8 
and z — A. In one dimension (d=l-fl), the roughness exponent, a = 1.5, exceeds unity, implying that the large 
scale steady state morphology of the growing interface is not self-affine in d=lH-l. The issue of whether the MH 
universality is only a crossover phenomenon or a true universality class is still debated [ ^0| . There are two lattice 
models to simulate MH growth. 

Das Sarma-Tamborenea (DT) Model - In this model |2|] a particle, after being deposited on a randomly chosen 
site, relaxes only to a nearest kink site i.e. it seeks only to increase the number of neighbors. 

Kim-Das Sarma (KD) Model - Choose a site i randomly and add a particle on i-1, i or i-l-1 (in d=l-|-l), depending 
on which site has a larger curvature V^/i = [h{x -I- 1) -I- h{x — 1) — 2h{x)] than that at the site i ||3l[]. If there is more 
than one satisfying the larger curvature condition, we choose one among them randomly. If the curvature of the site 
i is one of the larger curvatures, we add the particle at site i. 

5) Molecular Beam Epitaxy (MBE) Universality - The most relevant universality class in the context of con- 
served epitaxial growth is the Molecular Beam Epitaxy (MBE) universality (also goes by the name Lai-Das Sarma- 
Villain (LDV) universality) |Q . This is described by the non-linear version of the MH surface diffusion equation 
with 

Ambe = -KV^h + XiV^iVhf. 

The scaling exponents are known from a one-loop DRG calculation which gives a — {5 — d)/3, (3 — {5 ~ d)/{7 + d) 
and z = {7 + d)/3. There are two discrete models to simulate the MBE universality class. 

Lai-Das Sarma (LD) Model - This model is similar to the DT model described above with the difference that if an 
atom falls in a kink site, it is allowed to break its bonds and jump either down or up to the nearest kink site with the 
smallest step height |3^ . 

Kim-Das Sarma (KD) Model - This is the generalization of the KD model |^l[] described above. The only difference 
is that the linear curvature is replaced by a non-linear curvature =[h(x-l) + h(x+l) -2h(x)] - ^ [^ix — 1) — h{x)]'^ 

+ [h{x + 1) - h{x)]^ y 
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IV. COMPUTER SIMULATIONS OF CYCLICAL GROWTH 



A. Introduction 

1. Simulation methods 

To test our scaling hypothesis for cychcal processes we performed numerical simulations in d=l+l using specific 
discrete growth models described above. The system size in the simulations was changed between 128 to 4096 lattice 
spacings. Periodic boundary condition is employed so that columns i and i+L (L=system size) are equivalent. A 
typical cycle consisted of a deposition of 5-20 layers (average number of particle deposited per site) and desorption 
of between 10% to 100% of the deposited amount. The maximum number of cycles n varied between 500 — 10000 
to reach the saturation. To obtain good statistics we took average over 50-5000 independent runs, depending on the 
pairs of primary processes and the system size. The smaller the system the larger was the number of runs to obtain 
good result. 

2. Time dependence of surface roughness 

In a single process of growth or desorption, the roughness scales with time (Eq. (2)). For cycling, we have 
theoretically shown that time is replaced by the new scaling variable: the number of cycles n. One might ask: how 
does the cyclical roughness change with actual time leading to a scaling in terms of nl Below we discuss how the 
time {t) dependence of the surface roughness of cyclical growths can be studied to understand the emergent scaling 
behavior in terms of the number of cycles n. 

Consider a linear growth (or erosion) process. The structure factor for this can be given by Eq. (11). From that 
expression of S{q,t), the surface width = ^/ L'^^^S{q,t) can be expressed as |Q 

W\t)^W^{t) + Wl,,{t), (30) 

Wfiat{t) would be the roughness for growth induced on a flat initial substrate, and Wo is the contribution due to 
the width of the rough substrate surface. Since the total width is the sum of a decreasing {Wo) and an increasing 
{Wfiat) part, Competition between the two terms can make it increase or decrease from the initial roughness VFo(O) 
for sometime (eventually the roughness will exceed VFo(O)). 

In cycling two primary processes act alternatingly. The roughness generated by one process is taken as the initial 
roughness PFo(O) for the second process. In Fig. 3 the height profile is shown for the cyclic growth with Random 
deposition and EW desorption. EW dissolution smoothes the very rough surface produced by random growth in 
one cycle. Also note that the roughness increases with the number of cycles n. In Fig. 4 we show the actual time 
dependence of the simulated cyclical growth composed of two linear primary processes belonging to the EW and MH 
universality classes. In one cycle, MH growth increases the surface width and then EW surface relaxation smoothes 
out the surface to lessen the roughness. The average behavior of surface roughness in terms of cycles increases and 
gives rise to scaling. The scaling exponents are determined by the relevant (in the RG sense) terms of the primary 
processes. In our example (Fig. 4), MH/EW cycles are dominated by the surface relaxation of EW model and gives 
rise to EW exponents. The change in roughness within a cycle can be thought of fluctuations to the average behavior, 
and becomes less important as the number of cycles is increased (Fig. 4). Similar scaling (in terms of n) continues to 
hold for cyclic growths when one or both the primary processes are non-linear. 

3. Extraction of scaling exponents 

The growth exponent /3 is extracted for different system sizes L. The roughness W vs n is plotted on a log-scale 
and the slope of the best fitted straight line yields the exponent (3. The value quoted is from the largest L (once the 
effective (5 became close to the asymptotic value). From Ws{L) = W{L, oo), the saturation width dependence on L, 
the roughness exponent a can be calculated. Simulation results for different system sizes are used to plot In Wg ~ InL^ 
which is fitted to a straight line whose slope measures the exponent a. In some cases we checked independently the 
value of a from the scale dependence of the height-difference correlation function (from a log- log plot of A(r, t) vs r 
and fit to Eq. (4)). 
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B. Simulation Results 



1. Linear primary processes 

For linear primary processes, we looked at the possible pairwise combinations of RD, EW and MH universality 
classes using the absorption/desorption algorithms described earlier in this paper. When EW and MH is combined 
with RD, wc obtain EW and MH exponents respectively, because those arc the processes which generate correlations 
on top of the random growth. Fig. 5 shows, for the RD/EW cyclical process, that the value of /3 (asymptotically) 
is independent of the relative duration of the two primary processes. MH/EW cycles produce an asymptotic cyclical 
scaling with EW exponents f3 = 0.258(5) (Fig. 6) and a = 0.52(3) (Fig. 7). We have also performed data collapse 
(Fig. 8) to establish the validity of our scaling hypothesis. Asymptotic EW exponents in MH/EW cycles confirms 
that the surface relaxation of EW is the dominating effect when paired with MH surface diffusion (or growth on kink 
sites). Our theoretical analysis also predicts that the EW scaling exponents will be imposed in the MH/EW cycles 
because EW has a smaller dynamic exponent {z = 2) compared to that {z = 4) of the MH universality class. 

2. One non-linear process with a linear One 

To simulate non-linear processes we used two lattice realizations BD and RSOS (with equivalent results) for the 
KPZ universality class and KD algorithm for the MBE universality. First we combine the KPZ universality with the 
EW or the MH universality in a cycle to see the effect of the non-linear KPZ term to the scaling exponents. The 
exponents obtained are (3 = 0.311(5), a = 0.51(1) for EW/BD. The exponent /? increased slowly with the system size 
and the effective (3 reached a value close to the asymptotic one only for the largest system size(i = 4096). For the 
reverse process BD/EW wc obtained and /? = 0.322(5), a ~ 0.50(1) (sec Fig. 9). These asymptotic exponents are 
consistent with the KPZ (3 = 1/3 and, of course, with a = 1/2, which is the common value of EW and KPZ. To look 
at primary processes with different values of a, a DT (ai = 1.5) deposition with a ballistic dcsorption (a2 = 0.5) 
were performed. We obtain the asymptotic values of the exponents for MH/KPZ: (3 = 0.311(15) and a = 0.48(2), 
both consistent with the KPZ values. In all these cyclic processes the KPZ non-linearity (V/i)^ remained relevant 
with respect to the linear terms (V^/i or V^/i) and retain its scaling exponents. To show the opposite behavior where 
the non-linear process is not the dominant one we simulated MBE/EW cycles. From our approximate RG procedure 
for the cyclic process, it is clear that the fourth-order non-linear perturbation V^(Vft.)^ is irrelevant with respect to 
the linear EW term V^/i. In our simulation, when we allow surface relaxation of the EW model only to the nearest 
neighbor we get an effective exponent {(3 ^ 0.31) different than EW (3 even for the largest system size (L = 4096) 
we used. The reason may be that the next nearest neighbors of a chosen site affect the curvature dependent growth 
process in the simulation of MBE. When we consider next nearest neighbors for surface relaxation in EW process in 
MBE/EW cycles we get a = 0.50(2) and (3 = 0.251(3) consistent with EW values. 

3. Two non-linear processes 

Finally, we tried cycles consisting of two non-linear primary processes belonging to the KPZ and MBE universality 
classes. Simulations of RSOS/RSOS (note that they are not time-reversed images of each other because of the non- 
linearity) gave surfaces with KPZ scaling for Ti 7^ T2. For Ti = T2, however, EW behavior was found. This follows 
from the non-linear KPZ terms in the primary processes having the same magnitude but opposite signs. Hence, they 
exactly cancel each other in the coarse grained growth equations yielding an EW one. When we combine KPZ (using 
BD and KK) with MBE (KD) we exp(;ct to get KPZ scaling for cycles due to the dominant KPZ non-linearity (in 
d=l-|-l, where perturbative RG is applicable). Though both the processes have the same (3 (in 1+1), in our simulation 
of cycles we observe slow increase of beta with the system size when the BD model was used to simulate KPZ growth. 
The value of effective beta for the largest system size (L = 8192) did not reach the asymptotic value (1/3). Also, 
some sort of crossover behavior is observed which varies depending on the ratio of Ti and T2. The BD model itself 
is not very efficient to reach the asymptotic regime of KPZ growth and probably that is reflected in the simulations 
of cyclical growth. When We used the RSOS model for KPZ part of the cycle, the asymptotic scaling exponents 
(a = 0.50(1) and (3 = 0.335(5) ) were obtained for relatively small system sizes. 

We also run the respective absorption/absorption (or growth/growth) cycles with identical roughness behavior. 
Fig. 10 shows the results of DT/BD growth/growth cycling simulation. The scaling exponents {j3 = 0.31(2) and 
a = 0.48(2)) again acquire the corresponding KPZ values as expected. 
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V. EXPERIMENTS ON CYCLIC GROWTH 



Experiments of cyclical growth were performed by metal electrodeposition /dissolution of silver. The scaling behavior 
of surface roughness has been studied for electrochemical deposition [g5| Q |^ and dissolution processes 
separately. In cycling, however, metal electrodeposition, followed by a partial dissolution occurs. The substrate is 
plated for a specific period of time and then the current is reversed and the metal is dissolved from the substrate. 
This phenomenon is inherent in rechargeable batteries. 

Multiple cycles were carried out on vapor-deposited silver substrates, ranging from 1 to 20 cycles. The plating solu- 
tion contained 0.092M AgBr (silver bromide), 0.23M {NHa)2S20z (ammonium thiosulfate) , and 0.17M {NH4)2SO^ 
(ammonium sulfite). Each cycle consisted of plating for 5 min followed by 2.5 min of electrodissolution with a current 
density of 0.8mA/cm^. Image and scaling analysis were done using an atomic force microscope (AFM). Roughness 
was measured after n full cycles and after the deposition part of the cycles ( i.e. after n + p cycles with p — 2/3). A 
logarithmic plot of saturation rms height versus number of cycles, shown in Fig. 11, resulted in a straight line over 
two decades, indicating that the roughness scales with the number of cycles n, with P = 0.52, where (3 is the growth 
exponent. The roughness for cycling and deposition alone is compared. The value of /? is 0.52 for the cycling processes 
and 0.62 for deposition alone [Q, suggesting that processes such as the erosion of rough areas and the fiUing in of 
surface recesses are occurring more during cycling than in straight deposition. However, the saturation rms height 
was larger in magnitude for cycling than for straight deposition. 

The results suggest that cycling causes smoothing of the surface recesses and valleys by the dissolution of large 
peaks during the reverse plating process. 

VI. CONCLUSIONS 

In this paper we have discussed cyclical growth processes using dynamic scaling in terms of the number of cycles n 
(replacing the time variable t of a simple growth) . Given two primary processes we have described how to derive the 
scaling exponents of the combined cyclical process. MC simulations of different cyclic processes were carried out using 
various simple pairs of primary growth processes. Results of numerical simulations are entirely consistent with our 
theoretical understanding. For linear processes, the primary process with smaller dynamic exponent z always imposes 
its own scaling exponents on the combined cyclical process. The application of the RG method to cyclical processes 
was presented to determine cyclical scaling exponents in the case of non-linear primary processes. In particular, 
cyclical scaling exponents are independent of the relative duration (p or 1 — p) of the primary processes and also of 
the duration T of the cycle. However, both p and T do affect the amplitude of the roughness and also the crossover 
length beyond which the asymptotic exponents show up. First experimental results on electrodeposition- dissolution 
cycles are in accordance with our scaling hypothesis. 

Cyclic growth phenomena are widespread in Nature: The sea shores are shaped by cyclical high and low tides 
associated with the lunar motion around the Earth, all floral growth is subjected to daily changes in light and 
seasonal variations (which also affect many geological processes), etc. Thus the cyclical scaling theory introduced here 
might be found helpful in future investigations of many natural growth phenomena. 

Cyclical processes are also ubiquitous in technological applications. For example corrosion processes are strongly 
affected by seasonal weather conditions. The most likely application of our theory is for rechargeable batteries. 
Indeed, one of the breakdown mechanisms is due to the metal grown on one electrode reaching the other one thereby 
causing a shortcircuit. To test batteries today one has to run them their full lifetime. Using cyclical scaling, however, 
accelerated testing will become feasible. Results from measurements on fewer cycles for a short amount of time could 
be extrapolated to predict the battery lifetime under its realistic working conditions. 

Finally application to medicine may also be envisioned. The first one which comes to mind is the treatment of 
malignant tumors. They are known to have rough surfaces which follow the kinetic growth scaling laws [ p^ . Under 
a chemotherapeutic (or radiation) treatment they shrink and, hopefully, disappear. These treatments are always 
cyclical, and thus the cyclical scaling approach should faithfully describe the tumor's recession. We hope that using 
this approach will help plan the schedule and the dosage of the treatments such as to make them more efficient while 
minimizing the side effects. 
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FIG. 1. The roughness W vs t in the EW desorption process {log-log plot). (Inset: saturation roughness In VKs vs InL). 
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FIG. 2. The roughness W vs t in the KPZ (RSOS) desorption process {log-log plot). (Inset: maximal roughness InWs vs 
InL). 
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FIG. 3. Height profile for RD/EW cyclical growth. L=512 and 40 particles/site of Random deposition (solid line) and 20 
particles/site of EW dissolution (broken line) are used in one cycle. 
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FIG. 4. InVF (roughness) ~ Int of the MH/EW cyclic process. DT process increases the roughness and EW smoothes the 
surface in a single cycle, but the average roughness increases and gives rise to scaling (in terms of cycles). 
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FIG. 5. The roughness W vs number of cycles n in the RD/EW cychcal growth with different ratios of deposition Ti and 
dissolution T2 {log-log plot). System size of L=1024 and deposition of 20000 particles (fixed) are used. 
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FIG. 6. In (roughness) vs Inn (number of cycles) of the MH/EW cyclical process for different system sizes L. 
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FIG. 7. Height-difference correlation function A(r) vs distance r for the MH/EW cyclic growth is plotted after diflferent 
number of cycles {log-log plot). Inset: \nWs (maximal roughness calculated from Fig. 6) vs InL. 
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FIG. 8. Data collapse for MH/EW cyclical growth (using the exponents found from the graphs in Fig.4 and Fig.5.) clearly 
shows scaling in terms of cycles n. 
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FIG. 9. InW^ (roughness) vs hi n (number of cycles) of the simulated KPZ/EW cyclic growth for different system sizes L. 
Inset: roughness exponent a is extracted from its maximal values Ws for different L. 
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FIG. 10. The roughness W vs n of the DT/BD (growth/gowth) cyclical process {log-log plot). Inset: saturation roughness 
In Ws vs In L. 
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FIG. 11. The roughness vs number of cycles n in the electrochemical cyclical growth of Silver {log-log plot). 
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